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We present a simple on-chip electronic thermometer with the potential to operate down to 1 mK. 
It is based on transport through a single normal-metal - superconductor tunnel junction with rapidly 
widening leads. The current through the junction is determined by the temperature of the normal 
electrode that is efficiently thermalized to the phonon bath, and it is virtually insensitive to the 
temperature of the superconductor, even when the latter is relatively far from equilibrium. We 
demonstrate here the operation of the device down to 7 mK and present a systematic thermal 
analysis. 


I. INTRODUCTION 

On-chip electronic thermometry is an important part 
of modern research and commercial applications of nano¬ 
technology, and it has been studied already for several 
decades; see Ref. [1] and references therein. Many of 
these thermometers are based on tunnel junctions or 
quantum dots [2-4]. Temperature sensors based on nor¬ 
mal (N) and superconducting (S) metal tunnel junctions 
are used in a wide range of experiments [5-7] and appli¬ 
cations [8, 9]. An example of such a device is a primary 
Coulomb blockade thermometer (CBT) that is based on 
normal-metal tunnel junctions with an insulator as a 
tunnel barrier (NIN) [10, 11], where the electronic tem¬ 
perature can be obtained by measuring the smearing of 
the single-electron blockade. One more example is a 
SNS thermometer [12], whose critical current I c depends 
strongly on the temperature. Primary electronic ther¬ 
mometry has also been successfully demonstrated down 
to 10 mK using shot noise of a tunnel junction (SNT) [13- 
15]. Nowadays, a standard dilution refrigerator reaches 
a temperature of 5-10 mK, with a record of 1.75 mK 
[16, 17]. Nevertheless, a thermometer that has a modest 
structure and a simple but accurate temperature read¬ 
ing at sub-10-mK temperatures and does not require a 
complicated experimental setup is still missing. For this 
purpose, we present a normal-metal - insulator - super¬ 
conductor (NIS) junction that is widely used both as a 
refrigerating element and a probe of the local electronic 
temperature in different experiments and applications [5- 
9, 18-20]. The possibility to use the NIS junction at sub- 
10-rnK temperatures makes this thermometer suited for 
cryogenic applications at low temperatures. For instance, 
quantum information is a highly focused and rapidly de¬ 
veloping held in modern physics. For many realizations, 
such as superconducting and quantum-dot qubits, one 
needs to define a set of quantum states at low tempera¬ 
ture, that are well separated and well controlled and in¬ 


sensitive to noise and decoherence effects [21]. Several ex¬ 
perimental realizations of two-level systems [22-24] sug¬ 
gest that decreasing temperature further will increase the 
coherence times as well as improve charge sensitivity. We 
think that our thermometer will be interesting for a com¬ 
munity who is willing to discover new physics as well as 
improve already existing devices that require low tem¬ 
peratures for their proper functioning. The NIS ther¬ 
mometer is easy to operate compared to SNT [13, 14], 
and its thermalization is quite straightforward compared 
to CBT [11] due to the single-junction configuration and 
can be combined on chip with other solid-state devices. 
A measurement of the NIS current-voltage (/ — V) char¬ 
acteristic yields a primary temperature reading. 

In this paper, we study both experimentally and the¬ 
oretically an on-chip electronic thermometer based on a 
single NIS tunnel junction at sub-10-mK temperatures. 
We demonstrate the operation of the NIS thermometer 
down to 7.3 mK. In addition, we develop a thermal model 
that explains our measurement data and shows that self¬ 
heating effects remain negligible for temperatures down 
to 1 mK. 


II. THEORETICAL BACKGROUND 

Transport through a NIS junction has strong bias and 
temperature dependence. Near zero bias voltage, the cur¬ 
rent is suppressed due to the superconducting gap, A 
[25]. When biased at voltage V, current depends on the 
temperature due to the broadening of the Fermi distri¬ 
bution / n(E ) = (e £: / fe B T N _|_ i)- 1 j n normal metal 
with the temperature T)v and the Boltzmann constant 
Ub- The current can be expressed as [25] 

I = f + dEn s (E)[f N (E - eV) — f N (E + eV)}, 

2eH T 

(1) 
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where Rt is the tunneling resistance of the junction, and 
E is the energy relative to the chemical potential. 

In the superconductor, the Bardeen-Cooper-Schrieffer 
density of states is smeared and typically described 
by the Dynes parameter 7 expressed as n s {E) = 
15 le(u/Vu 2 — 1)| (see the Supplemental Material in Ref. 
[26]), where u = E/A(Ts) + * 7 , and Tg is the temper¬ 
ature of the superconductor. Possible origins of 7 in¬ 
clude broadening of the quasiparticle energy levels due 
to finite lifetime [27], Andreev current [28, 29], as well as 
photon-assisted tunneling caused by high-frequency noise 
and black-body radiation [26]. The typical experimental 
range of 7 for Al-based [30] tunnel junctions is 10 -4 to 
10 ~ 5 for a single NIS junction [26, 31], getting as low as 
10 -7 in SINIS single-electron transistors with multistage 
shielding [32]. 




Figure 1. In panel (a), we show the relative deviation of the 
present thermometer reading using method B with a numer¬ 
ically calculated I — V, Eq. (1). Sets of curves present the 
values of t for three 7 parameters: 10“ 7 , 2 . 2 X 10 -5 (the ac¬ 
tual value in the experiment), and ICE 4 shown as the red, 
blue, and green curves from left to right. For each 7 param¬ 
eter, temperatures 1, 3, and 7 mK are shown as dash-dotted, 
dashed, and solid lines, respectively. All curves are calculated 
using the parameters of the measured device with A = 200 
peV and Rt = 7.7 kD. The main panel (b) shows the mea¬ 
sured I—V characteristic (blue dots) together with full fit 
(solid red line) enlarged in the superconducting gap region. 
(Inset) Measured and calculated I — V curve on a wider volt¬ 
age scale at approximately 10 mK. 


One can determine Ty from a measured / — V curve 
using Eq. (1). As we show below, the self-heating 
of both N and S electrodes has a small effect on 
the / — V characteristic; thus, for now, we neglect 
these effects. Therefore, we assume temperatures to be 
small, ksTisr^s "C eV, A(Ts), and the superconduct¬ 
ing gap to be constant and equal to its zero-temperature 
value A. In this case, for eV < A, one can approximate 
Eq. (1) by 


I ~ / 0 e-( A ~°v)/k B T„ + 


'yV 

Rt^I - (eV/AF 


( 2 ) 


where Iq = \J2-kA ksTj^/2eRr [ 6 , 28]. Here the second 
term stands for the corrections to the I—V characteristic 


due to smearing. It leads to the saturation of the expo¬ 
nential increase of the current at low bias values. In the 
regime of moderate bias voltages, one can neglect this 
term and invert Eq. (2) into 


y= A, + Mv (3) 

e e 

This equation provides a way to obtain the electronic 
temperature by only the fundamental constants and 
by the slope of the measured I — V characteristic on a 
semilogarithmic scale as 


Tg{V) 


e dV 
~ks d(ln I) ’ 


( 4 ) 


This Eq. (4) allows us to use the NIS junction as a pri¬ 
mary thermometer, however, with some limitations. One 
can include the effects of 7 into Eq. (4) by subtracting the 
last term in Eq. (2) from the current I and obtain better 
accuracy. We do not take this approach here, since the 
main advantage of Eq. (4) is its simplicity as a primary 
thermometer without any fitting parameters. 

Next, we will compare the two methods used to ex¬ 
tract the electronic temperature from the measured I—V 
curves. In method A , we employ Eq. (1) and perform a 
nonlinear least-squares fit of a full I — V characteristic 
with Tjv as the only free parameter. The value of Tjy ob¬ 
tained in this manner, named T 4 , is not sensitive to 7 . 
Method B is based on the local slope of the I — V [see 
Eq. (4)]. The smearing parameter 7 has an influence on 
the slope of the I—V characteristic and, thus, induces 
errors in the temperature determination. The tempera¬ 
ture Ty is extracted as the slope of measured V vs In I 
over a fixed / range for all temperatures where Eq. (4) is 
valid. In the experiment, it is usually difficult to deter¬ 
mine the environment parameters precisely, but one can 
determine 7 from the ratio of Rt and the measured zero 
bias resistance of the junction. The I — V, which takes 
the 7 parameter into account [see Eq. (2)], gives indistin¬ 
guishable results from the ones obtained by method A. 

We evaluate the influence of the 7 parameter on the 
relative deviations of the present thermometer based on 
method B numerically, as shown in Fig. 1 (a). This devi¬ 
ation t is defined as the relative error t = (T^/Tjy) — 1. 
We show the values of t vs ln(ei?T^/A) for two extreme 
cases 7 = ICE 7 , ICE 4 , and for 7 = 2.2xlCE 5 extracted 
from the present experiment at temperatures of 1, 3, 
and 7 mK. The lowest bath temperature is 3 mK, and 
7 mK is the saturation of the electronic temperature in 
the current experiment. The larger the values of 7 and 
the lower the temperature, the higher the relative devia¬ 
tions. In addition, the range of the slope used to extract 
Tb shrinks with increasing 7 and with decreasing tem¬ 
perature [see, e.g., red curves in Fig. 1 (a)]. Thus, reduc¬ 
ing the leakage will significantly improve the accuracy of 
the device, especially towards lower temperatures. Pos- 
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sible avenues for suppressing 7 include improved shield¬ 
ing [32, 33] and encapsulating the device between ground 
planes intended to reduce the influence of the electromag¬ 
netic environment [26]. Finally, higher tunneling resis¬ 
tance of the junction decreases the Andreev current [34], 
We note that one can also use dV/d(lng) as a primary 
thermometer, where g = dl/dV is the differential con¬ 
ductance - typically a more precise measurement since it 
is done with a lock-in technique. Compared to Eq. (4), 
this method has the minimal deviation t reduced by at 
least a factor of 3.5 for T)v > 1 mK (6 for T)v > 10 mK), 
though exhibiting qualitatively similar dependencies on 
7 and T/v. 

III. EXPERIMENTAL REALIZATION AND 
MEASUREMENT TECHNIQUES 

Next, we describe the realization of the NIS thermome¬ 
ter that is shown together with a schematic of the ex¬ 
perimental setup in the scanning-electron micrograph in 
Fig. 2. The device is made by electron-beam lithography 
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Figure 2. A scanning-electron micrograph of the NIS device 
together with a schematic of the experimental setup. In the 
main panel, the S and N leads of the junction are visible, and 
underneath the pads, the ground plane of a square shape is 
indicated by a dashed line. Enlarged inset shows the actual 
tunnel junction where the S and N leads are shown in blue 
and brown, respectively. 

and the two-angle shadow evaporation technique [35]. 
The ground plane under the junction is made out of 50 
nm of Au. To electrically isolate the ground plane from 
the junction, we cover the Au layer with 100 nm of AlO^ 
using atomic-layer deposition. Next, we deposit a layer 
of ds = 40 nm of A1 that is thermally oxidized in situ. 
The last layer is formed immediately after the oxidation 
process by deposition of dw = 150 nm of Cu, thus, cre¬ 
ating a NIS tunnel junction with an area A = 380 nm 
x 400 nm. The geometry of the junction is chosen such 
that the leads immediately open up at an angle of 90° 
and create large pads with an area Am = A$ = 1.25 mm 2 
providing good thermalization. The S lead is covered by 
a thick nornral-metal shadow as shown in brown in the 
inset of Fig. 2, where the N and S layers are interfaced 
by the same insulating layer of AlO^, as the junction. 

The experiment is performed in a dilution refrigera¬ 
tor (base temperature 9 mK) where each of the sam¬ 


ple wires is cooled by its own separate Cu nuclear re¬ 
frigerator (NR) [36], here providing bath temperatures 
Tbath down to 3 mK. Nuclear refrigerator temperatures 
after demagnetization are highly reproducible and ob¬ 
tained from the precooling temperatures and previously 
determined efficiencies [11]. Temperatures above approx¬ 
imately 9 mK are measured with a cerium magnesium ni¬ 
trate thermometer which is calibrated against a standard 
superconducting fixed-point device. Since the sample is 
sensitive to the stray magnetic field of that applied on 
the nuclear refrigerator, this field is compensated down 
to below 1 G using a separate solenoid. The I— V curves 
(see Fig. 2 for the electrical circuit) are measured using a 
home-built current preamplifier with input offset-volt age 
stabilization [37] to minimize distortions in the I — V 
curves. 

Filtering, radiation shielding, and thermalization are 
crucial for obtaining a low 7 and low device temperatures. 
Each sample wire goes through 1.6 m of thermocoax, 
followed by a silver epoxy microwave filter [38], a 30-kHz 
low-pass filter, and a sintered silver heat exchanger in the 
mixing chamber before passing the A1 heat switch and 
entering the Cu nuclear stage. The setup is described in 
detail in Ref. [11] and is further improved here (see the 
Appendix for more details). 

IV. RESULTS AND DISCUSSION 

In Fig. 1 (b), the measured I-V characteristic in 
the superconducting gap region is shown by blue dots. 
The solid red line corresponds to the full fit based on 
method A. In the inset, we present the I-V character¬ 
istic at a larger voltage scale, used to extract Rt■ In 
Fig. 3 (a), the measured I-V characteristics of the NIS 
junction are shown in logarithmic scale by blue dots at 
Tbath — 100,...,3 mK from left to right. The full fits 
are shown as dashed red lines. The tunneling resistance 
Rt = 7.7 kfl and the Dynes parameter 7 = 2.2xl0 -5 
used in all these fits are determined based on the I-V 
characteristics shown in Fig. 1 (b) at high and low volt¬ 
ages, respectively. For the lowest temperatures, Tm from 
the nonlinear fit depends strongly on the superconduct¬ 
ing gap [39], making it difficult to determine the gap 
with high enough accuracy [40]. However, Eq. (1) gives 
a possibility to perform a nonlinear least-squares fit, and 
Eq. (3) gives a linear fit, where the parameters A and Tm 
are responsible for the offset and the slope, respectively. 
Therefore, at high temperatures (approximately 100 mK 
in the present experiment), one can narrow down the un¬ 
certainty in A such that Tm becomes essentially an inde¬ 
pendent parameter for the fits. Thus, the gap extracted 
from the high-temperature data using Eq. (1) remains 
fixed, A = 200 ± 0.5 /xeV, for all temperatures below 
100 mK. In addition, we show as solid black lines an ex¬ 
ponential / — V dependence corresponding to method B 
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Figure 3. Panel (a) shows measured I-V characteristics (blue 
dots), when Tbath is lowered from left to right together with 
fits as solid black and dashed red lines (see text). (Inset) 
Close-up of two I-V characteristics for temperatures of 10 
and 7 mK. The electronic temperatures extracted from both 
the full fit (red squares) of the I-V characteristics and their 
slopes (black triangles) are shown in (b). 


with a fitting range between 5 pA and 400 pA. The en¬ 
larged inset shows the I— V curves at temperatures of 10 
and 7 mK. The I-V characteristics presented in Fig. 3 (a) 
agree well with the theoretical expressions in Eqs. (1) 
and (2) (the latter form is not shown). In Fig. 3 (b), we 
show the electronic temperatures obtained using method 
A (red squares) and method B (black triangles) vs Tbath- 
Method A ( B ) shows a relative error in the electronic 
temperature up to 6 % (11 %). The error in method B 
is larger, as we neglect the influence of the 7 parameter. 

The lowest temperature obtained from the full fit is 
T$ = 7.3 mK with statistical uncertainty of 5 % at 
Tbath = 3 mK. The NIS temperature decreases slowly 
over time, arriving at 7.3 mK several weeks after the 
cooldown from room temperature. This suggests that 
internal relaxation causing a time-dependent heat leak, 
e.g., in the silver epoxy sample holder, is limiting the 


minimum temperature. Future improvements will em¬ 
ploy low-heat-release materials better suited for ultralow 
temperatures such as sapphire or pure annealed metals, 
e.g., for the socket and chip carrier, minimizing organic 
noncrystalline substances such as epoxies. 

V. THERMAL MODEL 

The total power dissipated in the device is equal to 
IV = Qnis + QniS’ w here Q% IS and Q% IS are the heat 
powers to the normal metal and to the superconductor, 
respectively. The heat released to the superconductor is 
given by 

Qih s = ^ J E s n s (E)[f N (E - eV) - fs(E)}dE, 

(5) 

where Es = E is the quasiparticle energy. To evaluate 
Qnisi one has to substitute Es by E n = (eV — E) in 
Eq. (5). Almost all of the heat is delivered to the su¬ 
perconductor in the measured (subgap) bias range, thus, 

Qnis ~ IV an d Qnis ^ Qnis- 

So far, we neglected all self-heating effects both in the 
normal metal and in the superconductor. To justify the 
nonself-heating assumption, we check numerically and 
analytically these self-heating effects. We sketch the an¬ 
alytical arguments in the Appendix. Here we state the 
main results obtained from the thermal model. 

Self-heating of the superconductor can take place due 
to the exponential suppression of thermal conductivity 
and the weak electron-phonon (e — ph) coupling, espe¬ 
cially at low temperatures. We find that the supercon¬ 
ductor temperature Tg stays below 250 mK in the subgap 
bias range | Vj < ^ and does not influence the thermome¬ 
ter reading. In this bias range and at Tbath = 3 mK, we 
estimate based on the numerical calculations the temper¬ 
ature of the superconductor Tg = 145 mK and the power 
injected to the superconductor is IV ~ 90 fW. At the 
same time, we evaluate the relative change of the slope 
to be small |t| < 5 x 10 -3 at I < 1 nA. In conclusion, 
the temperatures obtained from both methods A and B 
are affected by less than 0.5 % by self-heating of the su¬ 
perconductor [41]. In addition, the normal metal might 
get self-heated due to weak electron-phonon coupling and 
backflow of heat from the superconductor [42]. The in¬ 
fluence of the self-heating of the normal metal down to 
1 mK temperature affects the temperature obtained from 
both methods A and B by less than 0.5 % as well, as in 
the case of self-heating of the superconductor. 

VI. CONCLUSIONS 

In conclusion, we demonstrate experimentally the op¬ 
eration of an electronic thermometer based on a single 
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NIS tunnel junction. The thermometer agrees well with 
the refrigerator thermometer down to about 10 rnK and 
reaches a lowest temperature of 7.3 mK at Ttath = 3 mK, 
currently limited by a time-dependent heat leak to the 
sample stage. We discuss several possible improvements 
of the present device and experimental setup. Finally, 
we show that self-lreating in the normal metal and in the 
superconductor on the full I-V or its slope is negligible, 
paving the way for NIS thermometry down to 1 mK if 
the experimental challenges can be overcome. 
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APPENDIX 

A. Experimental techniques 

The setup described in Ref. [11] is improved as fol¬ 
lows. First, a ceramic chip carrier is replaced by silver 
epoxy parts which remain metallic to the lowest tem¬ 
peratures, allowing more efficient cooling. Further, the 
sample - previously mounted openly inside the cold-plate 
radiation shield together with the nuclear stages - is en¬ 
closed in an additional silver shield, sealed with silver 
paint against the silver epoxy socket, and thermalized to 
one of the Cu refrigerators (see Fig. 4). Finally, each wire 
is fed into the sample shield through an additional silver 
epoxy microwave filter. While previously saturating at 
10 mK or above [11], metallic CBTs have given tem¬ 
peratures around 7 mK after the improvements [38, 43], 
comparable to the NIS temperatures presented here. 


B. Estimates of the subgap conductance 

The Dynes parameter 7 discussed in the main text can 
be attributed to the higher-order processes such as An¬ 
dreev tunneling events. Assuming ballistic transport and 
an effective area of conduction channel A c h = 30 nm 2 
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Figure 4. Scheme of a dilution unit together with a NR. Ra¬ 
diation shields (not shown) are attached to the still and cold 
plate (approximately 50 mK). The resistor-capacitor (RC) fil¬ 
ters are 1.6 kfi/2.2 nF and 2.4 kfl/470 pF. The 21 NR plates 
are 0.25 x 3.2 x 9.0 cm 3 each, amounting to 64-g Cu per plate. 
The NRs cool as low as 0.2 mK. In the present experiment, 
the lowest Tbath used was 3 mK. Compared to Ref. [11], 
the improved setup depicted here features a Ag epoxy socket, 
a Ag epoxy chip carrier, and a second filtering stage with 
radiation-tight feedthroughs into an additional sample radia¬ 
tion shield. The abbreviations BW, B c and GND presented 
in the schematic stand for bandwidth, critical magnetic field 
and electrical ground, respectively. 


[44, 45], a simple estimate of subgap Andreev conduc- 
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tance reads oar = Rk/(8NRt ) = 8.5xl0 -5 in units of 
R^ 1 , where Rk is the quantum resistance, N is the effec¬ 
tive number of the conduction channels, and N = A/A c h, 
where A is the area of the junction. Alternatively, the 
estimate based on diffusively enhanced Andreev conduc¬ 
tance yields of corresponding dimensionless conductance 
is 7.5 xlO -5 . These values are of the same order of mag¬ 
nitude as in our experiment (7 = 2.2xlO -5 ) and fall in 
the range of earlier experiments [29]. 


C. Theoretical estimates for the relative deviations 
of the present thermometer 

The theoretical deviations of dV/d(\rtI) 
at 7 = 2.2xlO -5 numerically calculated from I — V, 
Eq. (1), are rather large, particularly at low tempera¬ 
tures [approximately 30% at 1 rnK; see solid blue curves 
in Figs. 1 (a) and 5]. Measuring the differential con¬ 
ductance g = dl/dV rather than current I significantly 
reduces the predicted deviations t g = T^ ope ' 9 /Tn — 1, 
where T^ ope ’ 9 — d Q^ g ^ e/k B . The minimum of these 
deviations gets broader and potentially reduces mea¬ 
surement noise since it is a lock-in measurement - overall 
strengthening the method B. In Fig. 5, for comparison 



ln(eR T l/A) 

Figure 5. The theoretical deviations t and t g of the ther¬ 
mometer reading using method B based on I(V) and g(V), 
respectively. These deviations are shown for temperatures 1, 
3, and 7 mK as dash-dotted, dashed, and solid lines, respec¬ 
tively. Parameters used are 7 = i2.2x 10 -5 , A = 200 /xeV, and 
Rt = 7.7 kfl are as in the actual experiment. 

we show two sets of curves for t g (thick purple curves) 
and t (thin blue curves) from left to right. Sets are 
calculated based on the experimental parameters for 7 
= 2.2xl0~ 5 , A = 200 /xeV and Rt = 7.7 kfl. Each 
set corresponds to the temperatures 1, 3, and 7 mK 
and is shown as dash-dotted, dashed, and solid lines, 
respectively. Here, the t set is identical to the set with 
7 = 2.2 x 10 ~ 5 that is shown in Fig. 1 (a). 


D. Self-heating of the superconductor 

First, we consider the self-heating of the superconduc¬ 
tor. We study the heat transport in the present geome¬ 
try (see Fig. 2) by a diffusion equation assuming thermal 
quasiparticle energy distribution [46, 47] 

-V(«sVT s ) = u s , (6) 

where we set a boundary condition near the junction 
- n in nerKS^Ts\junct = QnIs/ A i where dinner IS the in¬ 
ner normal to the junction. Thermal conductivity in the 
superconductor is 

6 . A -a 

K S = -7 2 e«i 0 r S ^, 7 

7T^ k B lg 

where Lq is the Lorenz number, and ctai = 3 x 10 ' 
(flm ) -1 is the electrical conductivity of the A1 film in 
the normal state [47], and we take into account the Tg 
dependence of the gap at low temperatures, A(T]g)/A ~ 
1 — yj2Trk B Tg /Ae -A / fcsTs . The absorbed heat is given 
by us = + Qtmp ■ Here, the first term is the electron- 

phonon power qf_ ph ~ Saz(T| - T^) exp(-A/k B T s ) 
[48], where Y,ai = 3 x 10 s WK - 5 m -3 is the material- 
dependent electron-phonon coupling constant. The 
phonon temperature T p is assumed to be equal to Tbath- 
Because of weak electron-phonon coupling, nearly all the 
heat is released through the (unbiased) normal-metal 
shadow (see Fig. 6 ) that acts as a trap for quasiparticles, 
qtrap- Here the conductance of the trap per unit area 

Junction . 



Figure 6 . The thermal diagram of the NIS thermometer. 
Present schematic does not reflect the real thicknesses of the 
materials. In this thermal model, we assume the normal- 
metal shadow that acts as the trap to be at Tbath- 

is the same as for the tunnel junction ot = 1 /(RtA). 
Therefore, the heat removed per volume by this trap qtrap 
can be calculated using Eq. (5) at V = 0, Tjv = Tbath and 
substituting Rt by ds/oT- Thus, the temperature of the 
superconductor Tg can be found as the solution of Eq. ( 6 ) 
in 2D in polar coordinates using radial approximation for 
the sample geometry and can be written as [47] 

^/2nk B T s /A(T s )e- A ^ /kBTs = aQ s NIS . (8) 

Here, we indeed assume Q^is ~ -TV, and a = 
sjne 2 G/{dscrAi \/^k B TsA 3 (Ts)) is a coefficient that de¬ 
pends on Tg, and dimensionless parameter G = ln *A / ' r ’ ( b ~ 
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2,..., 3 is logarithmically dependent on the sample ge¬ 
ometry [47]. Here, A is the relaxation length of the order 
of approximately 10, ..100 ^m, and ro = 2A/(7rci/v) ~ 
500 nm is the radius of the contact in the present de¬ 
vice. After substitution of all the parameters, we find 
that the superconductor temperature T$ does not influ¬ 
ence the thermometer reading staying below 250 mK in 
the subgap bias range |V| < ^. We estimate Tg to be 
approximately 145 mK in this bias range at T^ at i , = 3 mK 
corresponding to the power injected to the superconduc¬ 
tor as IV ~ 90 fW, and the quasiparticle density [47] 
as n qp = 0.3 pm~ 3 . In addition, we evaluate the rela¬ 
tive change of the slope to be small |f| < 5 x 10 -3 at 
/ < 1 nA. In conclusion, the temperatures obtained from 
both methods A and B are affected by less than 0.5 % 
by self-heating of the superconductor. 


area is A W i re = tt r 2 iT , e with a radius r W i re = 16 pm. 
The temperature obtained from both methods A and B 
is affected by less than 0.5 % by self-heating of the nor¬ 
mal metal down to a temperature of 1 mK. In addition, 
we can evaluate at low temperatures Tjv < Ts A /k B 
the maximum cooling at optimum bias voltage V op t ~ (A 
- 0.66 k B T N )/e [1] 


Q^IsiVopt ) 


A 2 
e 2 Rt 


-0.59 


2irk B Ts 


A 


x exp 


k B T N 

A 

k B T s 


3/2 


( 10 ) 


to be 90 pW at Tb at h = 1 mK. 


E. Self-heating of the normal metal 

For self-heating of the normal metal, one can solve the 
diffusion equation Eq. (6) taking into account the same 
boundary condition as above with all S indices replaced 
by N , where kn = LqcfcuTn is the thermal conductivity 
of the normal metal, and the electrical conductivity of Cu 
is assumed to be ac u = 5 x 10 7 (flm) -1 [49]. The heat 
absorbed in the normal metal is ujv = q^L p h + Qwire ■ The 
heat conduction through the gold bonding wires q£f ire is 
taken into account only at the point where it is attached 
to the normal-metal pad, whereas the electron-phonon 
interaction q£_ ph is effective in the full volume of the 
normal metal. The volumetric electron-phonon power is 
q»_ ph = £ Cu (T 3 - T 3 ), where £ c « = 2 x 10 9 W K~ 3 
m -3 is the electron-phonon coupling constant of copper. 
Here we consider the effect of the heat removed by the 
bonding wires on the temperature only in the normal 
metal, thus, q” ire = L 0 a Au (T% - Tg ath )/2L wire d N . The 
length of the gold bonding wire is L wire ~ 5 mm and a Au 
= 1.8 x 10° (flm) -1 is the electrical conductivity of gold 
measured at low temperatures. The thermal relaxation 
length in the normal metal [1] is 

In = (T^ 1 )- 1 y/(a c M(2Z Cu ). (9) 

We substitute p = 5 and Tb at h = 10 mK and obtain In = 
17.5 mm. Since all the dimensions of the present device 
are smaller than 1.5 mm, there is only a weak temper¬ 
ature gradient over the normal-metal electrode due to 
its good heat conduction, and the electron-phonon cou¬ 
pling weakens at low temperatures. By solving the heat- 
balance equation Qnis = Qe-ph+Qwire an d assuming no 
external heat leaks, one can calculate T)v. Here, the heat 
released through e — ph coupling is Q p _ ph = ^Nde-ph, 
where fijv = A n cIn is the volume of the N electrode. 
The heat released through N W i re = 2 bonding wires is 
Qwire = Qwire N wireA w i re d N , where its cross-sectional 
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